Interacting lattice electrons with disorder in two dimensions: Numerical evidence for 
a metal-insulator transition with a universal critical conductivity 
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The dc-conductivity of electrons on a square lattice interacting with a local repulsion in the 
presence of disorder is computed by means of quantum Monte Carlo simulations. We provide 
evidence for the existence of a transition from an Anderson insulator to a correlated disordered metal 
with a universal value of the critical dc-conductivity (Jdc,crit = (1.18 ± 0.06) /h at the transition. 
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The Coulomb interaction between the electrons and 
the presence of disorder both strongly affect the proper- 
ties of solids Namely, electronic correlations and 
randomness are separately driving forces behind metal- 
insulator transitions (MITs) due to the localization and 
delocalization of particles. While the electronic repulsion 
may lead to a Mott-Hubbard MIT the scattering of 
non-interacting particles from randomly distributed im- 
purities can cause Anderson localization 0, Q . The si- 
multaneous presence of disorder and interactions lead to 
further subtle many-body effects which raise fundamen- 
tal questions in theory and experiment not only in solid 
state physics [iHl, Is Hl^. b ut also in the field of cold 
atoms in optical lattices [13| . 

Accor ding to the scaling theory of Anderson localiza- 
tion 

0, [la 

non-interacting electrons in two spatial di- 
mensions [d = 2) are localized in the presence of disor- 
der. Hence, in the thermodynamic limit at zero temper- 
ature, there is no metallic state in d = 2. By contrast, 
the experimental observation of a metal-insulator transi- 
tion (MIT) in resistivity measurements on various high 
mobility heterostructure samples and Si-MOSFETs [Ig 
clearly indicates that interactions can turn an Ander- 
son insulator into a metal. Near-perfect scaling of the 
resistivity data [l3| was taken as evidence for the pres- 
ence of a quantum critical point between the metallic and 
the Anderson localized state [1, . Recent investigations 
based on a non-linear sigma model (NLcrM) for 
interacting electrons with disorder in the continuum con- 
firm the existence of a such quantum critical point which 
is characterized by a universal value of the dc-resistivity. 
Universal critical conductivities were also discussed in 
other two-dimensional systems, e.g., in connection with 
the transition from a superconductor (superfluid) to an 
insulator 13, 21 1, in the integer quantum Hall effect (i^ . 
and in graphene [l^. 

Numerical investigations of the interplay between dis- 
order and interactions usually address electrons on a lat- 
tice rather than in the continuum. Various approaches in- 
clude Hartree-Fock investigations in c? = 3 [2J| and c? = 2 



[25|, quantum Monte Carlo (QMC) simulations [26|-l28|, 
and dynamical mean-field theory [291432 1. In their QMC 
studies of two-dimensional lattice electrons Dcnteneer et 
al. |26l . |27[ indeed found a phase transition between an 



Anderson insulator and a metallic phase in accordance 
with experiment [5i] . There has also been the proposal of 
the MIT percolation transition [33[ . 

In this Letter we provide evidence through extensive 
QMC simulations that in the Anderson-Hubbard model 
in (i = 2 there exists a transition between a metallic phase 
and an Anderson insulator, and that this transition takes 
place at a value of the dc-conductivity crdc,crit which is 
essentially independent of the critical interaction, criti- 
cal disorder, and particle density. The computation of 
such a universal value of the critical dc-conductivity pro- 
vides an explicit link to results obtained from effective 
theories in the continuum 181. Indeed, numerical inves- 



tigations of microscopic lattice models can provide details 
of the properties of a system at a quantum critical point 
which are not accessible within effective perturbative ap- 
proaches. 

Our investigation of interacting electrons in the pres- 
ence of disorder is based on the Anderson-Hubbard 
Hamiltonian on a square lattice 



H = T{e,} 



(1) 



Here 



<lj><T 



^{ei- fi)n,„, (2) 



is the single-electron part where c]^, (cio-) are fermion 
creation (annihilation) operators for site Ri and spin cr, 
i^ia = c[a'^i<y is the Operator for the local density, ^ de- 
notes the chemical potential, and t is the hopping am- 
plitude for electrons between nearest neighbor sites. The 
local energies are random variables which are sam- 
pled uniformly from the interval [—A/2, A/2]; hence the 
width A characterizes the strength of the disorder. The 
interaction is assumed to be repulsive {U > 0). The 
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model is solved numerically using determinantal QMC 
(DQMC) dl where the interval [0,^] {p = l/kBT) is 
partitioned according to (3 ~ L x At, with At as the 
size of a small step in the imaginary time direction, and 
L as the number of imaginary time slices. The parti- 
tion function Z is then decomposed according to the 
Suzuki- Trotter formula [35[ . In the next step, a Hubbard- 
Stratonovich transformation is performed whereby the 
interaction problem is reduced to non-interacting elec- 
trons in the presence of infinitely many fluctuating fields 
described by Isin g v ariables on every space- (imaginary) 
time lattice site [36[. The electrons can then be inte- 
grated out. The calculation of quantities such as the 
Green function, electronic density and two-particle cor- 
relation functions proceeds with Monte Carlo sampling 
of the various configurations of the Ising degrees of free- 
dom. The hopping integral t sets the unit of energy and 
the simulation now contains three independent energy- 
scales: the disorder strength A, the interaction strength 
U, and the temperature T. 

To evaluate the dc-conductivity, we compute the elec- 
tronic current density operator 



ieat 



E 



.a) , (3) 



where e^^ denotes a translation in x-direction by a lattice 
constant a. This leads to the time-dependent current 
density operator 



Jx 



(4) 



where t is the imaginary (Matsubara) time. The 
position-space Fourier transform of the current operator, 
ja;(q, t), is then used to calculate the current-current cor- 
relation function 



Aa;a;(q,T) (jx (Q, t) jj; (-q, T = 0)). 



(5) 



Within linear response theory, the dc-conductivity is ob- 
tained from 



CTdc = lini 



ImA^ 



0,w) 



(6) 



The current-current correlation function in Matsubara 
time is related to the imaginary part of the current- 
current correlation function in real frequency through the 
integral transform 



J^xxiq, t) 



dio 



n 1 



(7) 



DQMC simulations can compute A^j;(q, t), but to de- 
termine CTdc it is necessary to obtain ImAa;2;(q, w). For 
low enough temperatures, the exponential decay of the 
bosonic kernel K{lo,t\ /?) = j^l^./j^ for t = /3/2 ensures 
that the integral contributes only for small w, where the 
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FIG. 1: (Color online) Curves of the dc-conductivity vs. tem- 
perature at electron density n = 0.5 (quarter filling) and in- 
teraction i7 = 3 on a 10 X 10 square lattice computed for 
different values of the disorder strength A (see inset). 



substitution arising from linear response, eq. ©, is valid. 
Replacing t by /3/2 and ImAa;a;(q, w) by WCTdc, the inte- 
gral can be carried out analytically and yields the dc- 
conductivity (37I I as a function of temperature for differ- 
ent values of the interaction U and disorder strength A: 



0,T 



(8) 



In the following discussion, we set i = 1. 

The conductivity data is averaged over 10 disorder re- 
alizations at high temperatures (T = 1,0.5,0.333,0.25), 
up to 80 disorder realizations for intermediate temper- 
atures {T ~ 0.2,0.167), and up to 100 disorder realiza- 
tions for the two lowest temperatures (T = 0.125,0.1). 
In Fig. [1] the dc-conductivity is shown as a function 
of T for several values of the disorder strength A. Ini- 
tially, when the value of the disorder strength is less than 
about A = 10, the slope of the conductivity curve at low 
temperatures is negative (i.e., the conductivity decreases 
with increasing temperature), implying that the system 
is metallic. As the disorder strength is increased, the 
low temperature conductivity develops a positive slope, 
which is the signature of an insulator. Since the system is 
far from half-filling, such that a Mott- Hubbard MIT does 
not occur, these results indicate a transition between a 
metallic and an Anderson localized state. 

There are two sources of statistical error in this analy- 
sis: one due to the QMC simulations, the other due to the 
disorder averaging. For all parameter sets studied here, 
the intrinsic QMC error for any given disorder realiza- 
tion is much smaller than the error arising from different 
disorder realizations. Since the error bars in Figs. [1] [5] 
are of the order of, or smaller, than the symbols they arc 
not shown. 

On the basis of Fig. [T] neither the critical disorder 
strength nor the value of the dc-conductivity ddc at the 
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FIG. 2: (Color online) Plot of the critical conductivity a^c 
vs. disorder strength A for three temperatures and the same 
values of U and n as in Fig. [1] The well-defined crossing point 
determines the critical disorder strength as Ac = 9.8 . The 
conductivity at the critical disorder has the value crdc.crit = 
l.We'^/h. 



TABLE I: The six parameter sets of the interaction strength 
U and electron density n employed in our investigation of 
correlated electrons in the presence of disorder on a square 
lattice are listed together with the computed critical disor- 
der strength Ac and the critical dc-conductivity crdc,crit at 
the transition between a disordered metal and an Anderson 
insulator. 



Index 


U 


n 


Ac 


(Tdccrit (e'^/h) 


a 


1.0 


0.3 


6.8 


1.19 


b 


2.0 


0.3 


7.8 


1.07 


c 


3.0 


0.3 


8.6 


1.19 


d 


2.0 


0.5 


7.9 


1.26 


e 


3.0 


0.5 


9.8 


1.19 


f 


3.0 


0.6 


10.5 


1.19 



critical point can be determined accurately. In Fig. [2] 
we therefore plot iTdc as a function of the disorder A 
for the three lowest temperatures simulated here, i.e., 
T = 0.167, 0.125, 0.1. When A < Ac, the dc-conductivity 
increases with decreasing temperature (metallic behav- 
ior), while for A > Ac, the conductivity decreases 
with decreasing temperature (insulating behavior). The 
three curves shown in Fig. [2] display a well-defined cross- 
ing point, at which the dc-conductivity is independent 
of temperature, thereby marking the critical point for 
the MIT. From the location of the crossing point one 
can read off the value of the critical disorder strength 
Ac. For U = 3.0 at quarter-filling {n = 0.5) we find 
Ac = 9.8, while the value of the critical conductivity 
is CTdccrit = 1.19e^//i. We will use this technique to 
evaluate the critical disorder strength and the critical 
conductivity for five more parameter sets {U, n) listed in 
Tabic m The results collected in this table can be sum- 
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FIG. 3: (Color onhne) Graphical presentation of the spread 
of the parameters Usted in Table |T] in percentage relative 
to p, where p is one of the quantities U, n, Ac and 
o"dc,crit- Here p defines the arithmetic mean of the pa- 
rameter p across the six parameter sets a-f. For example, 
U = ^ (1.0 4- 2.0 + 3.0 + 2.0 + 3.0 + 3.0) = 2.33. The quan- 
tity U does not have a physical meaning, but is only a refer- 
ence point to display the relative spread of the U values used 
in the six parameter sets; the same holds for the parameters 
n and Ac. By contrast, o"dc,crit varies very little across the 
parameter sets and clusters around 1.18e^/h. Therefore the 
quantity crdc,crit can be accorded physical meaning. 



marized as follows. In spite of the strong variation of 
the microscopic input parameters n and the disorder 
strength Ac at the transition between the metallic and 
the Anderson- insulating state in c? = 2, the associated 
critical dc-conductivity, crdc,ciit, is found to be essentially 
independent of these input parameters. The results are 
presented graphically in Fig. [3l where the value of the 
critical conductivity is seen to cluster around the value 
CTdcxrit = (1.18 ± 0.06) e^//i. This provides evidence for 
the existence of a universal value of the critical conduc- 
tivity. Indeed, the NLuM, in which the dc-conductivity 
appears as a coupling constant, predicts a universal 13811 
critical value of the conductivity (Tdc.crit ^ 1.06e^//i jl8l |. 
in close correspondence with our result. In obtaining this 
estimate, we assumed that the number of valleys appro- 
priate for our work is n^, = 1 [39| . Thus our results estab- 
lish a link between the microscopic Anderson-Hubbard 
model, and the low-energy effective theory provided by 
the NLcrM for the metal-insulator phase transition in two 
dimensions. 

In summary, quantum Monte-Carlo simulations of in- 
teracting lattice electrons in the presence of disorder in 
d = 2 provide clear evidence for a transition from metal- 
lic to insulating behavior as the disorder strength is var- 
ied. At the transition the value of the dc-conductivity 
is found to be given by CTdc,crit = (1-18 ± 0.06) e^/h, im- 
plying that the critical dc-conductivity is essentially in- 
dependent of interaction strength, electron density and 
the critical disorder strength. This points towards the 
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existence of a universal critical dc-conductivity. We ob- 
tained qualitatively similar results from investigations 
where site-disorder is changed to bond-disorder, the de- 
tails of which will be published elsewhere. 
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